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The quantum efficiency in the transfer of an initial excitation in disordered finite networks, mod¬ 
eled by the fc-body embedded Gaussian ensembles of random matrices, is studied for bosons and 
fermions. The influence of the presence or absence of time-reversal symmetry and centrosymme- 
try/centrohermiticity are addressed. For bosons and fermions, the best efficiencies of the realizations 
of the ensemble are dramatically enhanced when centrosymmetry (centrohermiticity) is imposed. 

For few bosons distributed in two single-particle levels this permits perfect state transfer for almost 
all realizations when one-particle interactions are considered. For fermionic systems the enhance¬ 
ment is found to be maximal for cases when all but one single particle levels are occupied. 
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I. INTRODUCTION 

An important question on complex quantum systems, 
which remains open to a large extent, addresses the con¬ 
ditions to have robust efficient transport of excitations 
across a disordered finite network [l]-0|. Whereas it is 
well known how to define a Hamiltonian system where 
perfect transport is obtained Si , the fact that such 
system requires the precise specification of a large num¬ 
ber of parameters makes it difficult to achieve in practice. 
This is the sense of robustness above: statistical changes 
in the parameters should lead to small fluctuations that 
preserve good efficiencies, instead of large transmission 
fluctuations. Clearly, the number of control parameters 
should be as small as possible. 

Inspired by Ref. [b], we study here the distribution of 
the transport efficiencies of an initial localized excitation 
in a disordered network of l sites, which is modeled by 
random Hamiltonian that includes many-body interac¬ 
tions, considering both fermions and bosons; to the best 
of our knowledge, quantum efficiencies of this type of dis¬ 
ordered networks have not been considered. In general, 
this question is of interest in a variety of fields, includ¬ 
ing understanding photosynthetic light-harvesting com¬ 
plexes 0 such as the Fenna-Matthews-Olson (FMO) 
complex [8|, or in quantum communication protocols 
across quantum spin chains nm. In either case, a 
model with few-body interactions seems realistic. 

The transport efficiency from an input state |in) to an 
output state |out), is quantified as the maximum tran¬ 
sition probability achieved among theses states within a 
time interval [0, T], The transport efficiency is defined 
asi 

^in,out = max |(out|U(f)|in)| 2 . (1) 

[0,T] 
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Here, U{t) is the unitary quantum evolution associated 
with the Hamiltonian of the system and T is a reason¬ 
able time scale (h = 1). The system is said to have 
perfect state transfer (PST) when 'Pm, out = 10- The 
nodes of the network are the basis states of the Hilbert 
space where the initial excitation is localized. We de¬ 
fine the Hilbert space by distributing n spinless particles 
(bosons or fermions) in l single-particle states. The n- 
body Hamiltonian that we consider consists of a random 
fc-body interaction among the n-particles states, with 
1 < k < n; the non-zero Hamiltonian matrix elements are 
related to the links of the network. This matrix model 
is known as the (bosonic or fermionic) k- body; embedded 
Gaussian ensemble of random matrices [lj, |14j. 

The motivation for choosing this matrix model comes 
from the fact that this ensemble displays correlations 
among the matrix elements 1111,13 , in the sense that 

the number of independent random variables is usually 
smaller than the number of independent links of the net¬ 
work. In addition, the bosonic ensemble for the spe¬ 
cific case where the bosons are distributed in two single¬ 
particle levels displays a systematic appearance of dou¬ 
blets in the spectrum [17]. Thus, one of the design prin¬ 
ciples required in Ref. (6), the existence of a dominant 
doublet, is automatic fulfilled. The extension to fermions 
arises naturally, and is also motivated by the transport 
processes within the FMO system [?J. 

As in previous studies jl8|, centrosymmetry p~9| is con¬ 
sidered and we find that it increases dramatically the best 
efficiencies among certain pairs of distinct sites, in com¬ 
parison to a non-centrosymmetric interaction; we also 
study in this paper the case where time-reversal sym¬ 
metry is broken, and the corresponding case including 
centrohermiticity [20j . Here, as usual, the benchmark 
for good efficiencies is 95%. In particular, for k = 1 
and few bosons distributed in two single-particle states, 
we find that PST is obtained for a large fraction of the 
best efficiencies of the ensemble when centrosymmetry 
or centrohermiticity are included; for larger values of k 
the best efficiencies are reduced. In the case of fermions 
and including centrosymmetry or centrohermiticity, we 



2 


show that below half-filling (n/l < 1/2, where n is the 
number of fermions and l the number of single-particle 
states) the best efficiencies are obtained for fc = 1; above 
half-filling, our numerics indicate that the best efficien¬ 
cies correspond to k ~ 1/2 for n = l — 1. In both cases, 
the benchmark is achieved with non-zero probability. 

The paper is organized as follows: in Sect.[II]we define 
the matrix model considered and centrosymmetry, which 
is introduced at the one-particle level and then extended 
to the k- and n-particle spaces. In Sect. mu we describe 
the results for the distribution of the best efficiencies of 
each realization of the ensemble for bosons occupying 
two-single particle states, with or without the addition 
of centrosymmetry or centrohermiticity; Sect. IIVI ana¬ 
lyzes the fermionic case along the same lines. Finally, in 
Sect. IV] we present the conclusions of this work. 

II. THE EMBEDDED GAUSSIAN ENSEMBLES 
AND CENTROSYMMETRY / 

CENTROHERMITICITY FOR INTERACTING 
MANY-BODY SYSTEMS 

We introduce the fc-body embedded Gaussian ensem¬ 
bles of random matrices for bosons and fermions, for the 
cases of orthogonal (0 = 1) and unitary (0 = 2) symme¬ 
tries, following Ref. [h| ■ We consider a set of l degenerate 
single-particle states | j), with j = 1,2,...,/. The associ¬ 
ated creation and annihilation operators for fermions are 
Oj and a : j . and hj and bj for bosons, with j = 1 
These operators obey the usual (anti)commutation rela¬ 
tions that characterize the corresponding particles. We 
define the operators that create a normalized state with 
k < l fermions from the vacuum as ibl = ib\ ■ = 

J~[k =1 , with the convention that the indexes are or¬ 

dered increasingly ji < ji < ■ ■ ■ < jk (a simplifies the 
notation for these indexes); the corresponding annihila¬ 
tion operators are ipk-, a = (V4- a )^- Likewise, the fc-boson 

states are given by \l- a = \), . jk = N a n* = i b ] s > where 

again ji < ji < ■ - • < jk ■ Here, N a is a factor that guar¬ 

antees the normalization to unity of xj-. a |0): if the index 
j repeats kj times, N a contains a factor (fcj!) -1 / 2 . 

The random fc-body Hamiltonian for fermions reads 

= (2) 

a ,7 

a similar equation holds for bosons replacing V4- Q by 

Xk-a- I 11 ®q. © the coefficients v^ are random dis- 
trilDuted independent Gaussian variables with zero mean 
and constant variance 

W !fa, 7 W fefa', 7 ' = v 0 (^a,7 ,( ^a',7 + (3) 

Here, 0 is Dyson’s parameter that accounts for the pres¬ 
ence (0 = 1) or absence (0 = 2) of time-reversal sym¬ 
metry (2lj| . the bar denotes ensemble average, and we set 
Vo = 1 without loss of generality. 


The Hamiltonian H^ acts on a Hilbert space spanned 
by distributing n > fc particles on the l single-particle 
states. Then, a complete set of basis states is given by 
the set V/LalO) for fermions (with l > n), and xj,. a |0) for 
bosons. The dimension of the Hilbert spaces are, respec¬ 
tively, Np = and JVg = ( Z+ ™^ 1 )- This defines the 
fc-body embedded Gaussian ensembles of random matri¬ 
ces mm. 

By construction, the case k = n is identical to the 
canonical ensembles of Random Matrix Theory [ 2 J, i.e., 
to the Gaussian Orthogonal Ensemble (GOE) for 0 = 1 
or the Gaussian Unitary ensemble (GUE) for 0 = 2. For 
k < n, the matrix elements of Hjf' 1 may be identical to 
zero and display correlations. The former property ap¬ 
pears whenever there are no fc-body operators that link 
together the n-body states, e.g., for k < n/2. Correla¬ 
tions arise because matrix elements of Hjf' 1 not related 
by symmetry may be identical. One of the notorious 
consequences of this is, for the bosonic ensemble in the 
dense limit (n fc,L for fc and l fixed) that the ensem¬ 
ble is non-ergodic [l6j]- I n particular, for l = 2, fc n 
and 0 = 1, the spectrum displays a significant number of 
quasi-degenerate states US¬ 
As mentioned above, centrosymmetry is an important 
concept for an optimal efficiency [l^|- A symmetric NxN 
matrix A is centrosymmetric if [ A , J] = 0, where J is 
the exchange matrix Jij = 5i,N-i +1 [3; f° r complex 
hermitian matrices, centrohermiticity is defined when 
JA T J = A [13- 

Imposing centrosymmetry to the fc-body embedded en¬ 
sembles is subtle. It can be introduced either at the 
one-particle level, which is the core for the definition of 
the fc- and n-particle Hilbert spaces, or at the fc-body 
level, where the actual (random) parameters of the em¬ 
bedded ensembles are set, or at the n-body level, where 
the system evolves. Considering a more realistic descrip¬ 
tion which includes a one-body (mean-field) term and 
a two-body (residual) interaction, H = Hk -1 + Hk= 2 , 
it seems unnatural to define a specific transformation 
for each term separately. Hence, we shall introduce it 
in the one-particle space, and compute how is it trans¬ 
ferred to the fc-body and n-body space, which depends 
on 0 as well. Then, in the one-particle space we de¬ 
fine Ji\j) = \l — j + 1) for j = 1whose matrix 
representation in the one-body basis is precisely the ex¬ 
change matrix. For two fermions, we define j 2 = 

Jia' h J\a' h = -ipl. l _ j2+l l _j 1+1 - In the last equality we 
keep the convention that the indexes are arranged in in¬ 
creasing order; then, the fermionic anticommutation re¬ 
lations impose a global minus sign, which can be safely 
ignored. This is generalized for fc particles as 

k 

■'"‘■in . n = II J i4 = > ( 4 ) 

S = 1 

where we have dropped any global minus sign. We note 
that in general the matrix Jk, as defined by Eq. ®, is 
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not an exchange matrix, i.e., the matrix with ones in 
the counterdiagonal and zeros elsewhere. This follows 
from the possible existence of more than one state that 
is mapped by Jk onto itself; in this case, we shall say that 
Jk is a partial exchange matrix. As an example, consider 
the case of fermions distributed in l = 4 single-particle 
states; for k = 2 the fc-particle space has dimension 6 . 

Then, J 2 if)l - 2 ,3 = 2 - 2,3 and j ^ 2 -ia = ^ 4 ; i, 4 i ignoring 

the minus signs mentioned above. Then, the entries in 
the J 2 matrix elements for these basis states are 1 in the 
diagonal. In contrast, for the case l = 4 and k = 1,3 the 
resulting matrices J\ and J 3 are exchange matrices. 

For bosons Jk is defined as for fermions using Eq. a 
with the corresponding change of the creation operators. 
In this case, again, Jk may not be an exchange ma¬ 
trix. As an example consider k = 2 and l = 3; then, 
we have J 2 x 2 - 22 = N a (Jib \) 2 — X 2-2 2 > and likewise 

J 2 X 21 3 = X 20 3 ! which shows that there are more than 
one basis states mapped onto themselves under J 2 . Yet, 
for the special case l = 2 that we study below, Jk is 
an exchange matrix for all fc, and Hk (in the ?r-particle 
space) inherits the centrosymmetry (or centrohermitic- 
ity) of Vk, i.e., [ J n , Hk] = 0, with J n the exchange matrix 
of appropriate dimensions. 

III. STATISTICS OF THE TRANSPORT 
EFFICIENCY FOR BOSONS DISTRIBUTED IN 
l = 2 LEVELS 

We analyze here the transport efficiency for 

bosons distributed in l = 2 single-particle states, con¬ 
sidering an arbitrary initial state xlv.jj and an arbitrary 
final one Xn-v We address separately the cases where 
the fc-body random interaction is invariant or not with 
respect to time-reversal (/? = 1 or 0 = 2 ), and also when 
centrosymmetry (or centrohermiticity) is or not addition¬ 
ally imposed. In the simulations described below, we 
considered for concreteness the case with n = 9 bosons, 
the corresponding Hilbert space dimension is Nb = 10, 
the number of realizations of the ensemble is 2000 , and 
T = 15 in Eq. (JT]). 

We discuss first the case where centrosymmetry or cen¬ 
trohermiticity, respectively for 0=1 and 0 = 2 , is ab¬ 
sent. The distributions of the efficiencies V^, v of the 
ensemble for any combination of the input and output 
states is in general rather broad, with negligible contri¬ 
butions to efficiencies above 0.95. More interesting is to 
focus on the distribution of best efficiencies V of each 
realization of the ensemble. The results are presented in 
Fig. [I] The empty histograms correspond to 0 = 1 and 
the shaded ones to 0 = 2. With respect to time-reversal 
invariance, the case 0 = 2 seems to yield marginally bet¬ 
ter efficiencies for a given value of fc, in the sense that the 
mean value of the distribution attains larger values and 
the distribution is somewhat narrower. Yet, we notice 
that the largest values of V are slightly dominated by 
the 0 = 1 case, except for k = 1. Regarding the depen¬ 


dence on k , the case k = 1 notably distinguishes itself 
as a special one. First, the distribution of V is domi¬ 
nated by a peak at 0.4, is wider than for other values of 
k , and displays the largest number of realizations with 
efficiencies larger than the benchmark value 0.95 (verti¬ 
cal dashed line). Apart from this case, other values of 
k display rather poor efficiencies, with smaller probabil¬ 
ities of reaching the benchmark value if they do it at all. 
Actually, the distributions for intermediate values of k 
(4 < k < 6 ) seem to be marginally better. In general, the 
states Xn ;A t and xt;v that yield the best efficiencies dif¬ 
fer by one boson in the occupation of each single-particle 
state, except for k = 9 where the best efficiency is uni¬ 
formly distributed for all pairs of states. This statement 
is consistent with the fact that k = n corresponds to 
a GOE or GUE. Clearly, these unconstrained embedded 
Gaussian ensembles yield poor transfer efficiencies, ex¬ 
cept maybe for k = 1 . 

We turn now to the best efficiencies of each realiza¬ 
tion of the ensemble when centrosymmetry or centroher¬ 
miticity is imposed. As mentioned above, for the specific 
bosonic ensemble with 1 = 2 , imposing centrosymmetry 
or centrohermiticity at the one-particle level carries over 
to the k- and n-particle spaces; hence, [Hk , J n ] = 0, with 
J n the Nb x Nb exchange matrix. In Fig. [2] we present 
the best efficiencies V of each realization of the ensem¬ 
ble for this case. In comparison to the results presented 
in Fig. [I] centrosymmetry (centrohermiticity) enhances 
dramatically the best efficiencies. In particular, the case 
k = 1 is remarkable since it displays perfect state trans¬ 
fer for all realizations of the ensemble for 0 = 2, while 
for 0=1 there is a small remnant of a peak around 0.4. 
Quantitatively, for 0 = 1 about 92% of the realizations 
exhibit efficiencies that are larger than the benchmark 
value, while for 0 = 2 this number is well above 99%; the 
mean of the distributions is 0.9556 and 0.9978, respec¬ 
tively. Increasing k shifts the mean value of the distribu¬ 
tion towards smaller values and widens the distributions, 
being the effect larger for 0 = 2 than for 0 = 1. Whereas 
for k = 2 some realizations still display almost PST (last 
bin of the histograms), this decreases for increasing val¬ 
ues of k. For 0=1, all values of k exhibit realizations 
whose best efficiencies are above the benchmark value, 
with better efficiencies appearing for smaller values of k. 

The states xhyj. and Xn-,v that display the best efficien¬ 
cies are those linked by centrosymmetry, i.e., JnXnui = 
Xn „- While for k = n all centrosymmetry-related pairs 
participate uniformly, for other values of k there is a sig¬ 
nificant dominance of the pair of states where all bosons 
are initially in either of the single-particle states. From 
the perspective of the network, these states are located 
at the edges. 

As mentioned above and manifested in Fig. [2Ka), the 
case k = 1 is quite special. Indeed, for k = 1 there are 
only three (0 = 1 ) and four ( 0 = 2 ) independent ran¬ 
dom matrix elements; centrosymmetry (centrohermitic¬ 
ity) imposes that they become only two (0 = 1 ) and three 
(0 = 2). The consequence of this is that the random 
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Figure 1. Normalized distributions of the best efficiencies V of each realization of the bosonic embedded Gaussian ensemble 
with no centrosymmetry (or centrohermiticity) for l = 2, n = 9: (a) k = 1, (b) k = 2, (c) k = 3, (d) k = 4, (e) k = 5, (f) k — 6, 
(g) k = 7, (h) k = 8, and (i) k = 9. The empty histograms correspond to /? = 1 and the shaded (violet) to (3 = 2. The red 
vertical line indicates the 95% benchmark for the efficiency. Note that the scales for the k = 1 case are different from the rest. 
The inset in (a) is an enlargement of a region close to the benchmark value, showing the non-zero probability of having PST. 


many-body centrosymmetric/centrohermitian Hamilto¬ 
nian has a constant main diagonal, which is propor¬ 
tional to the total number of bosons, and a second- 
diagonal whose matrix elements are proportional to 
(rii + 1 ){Nb — rii — 1), where n» = 0,..., n is the bo¬ 
son occupation number of one of the single-particle lev¬ 


els for the many-body state in question. Except for the 
constant diagonal and the random weights involved, for 
f3 = 1 this matrix model has been shown to exhibit 
PST [22|. It can be shown that the eigenvectors of both 
models are identical, and the eigenvalues are related by 
a linear function. This ensures the occurrence of PST 
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Figure 2. Same as Fig. |T] where we impose centrosymmetry or centrohermiticity for the embedded Gaussian ensembles. Note 
that for k = 1 and /3 = 2, almost all realizations yield almost PST. 


as long as the off-diagonal matrix element of the random 
one-body interaction matrix is different from zero. Yet, 
the time scale where the PST occurs may be quite long, 
since it is inversely proportional to the absolute value 
of the off-diagonal matrix element of the random one- 
body matrix; this explains the occurrence of other values 
for the efficiency, as shown in Fig. [2fa). Our results ex¬ 
tend the validity of these statement to the broken time- 
reversal symmetric case and, more important, illustrate 


that they are robust under certain random centrosym- 
metric/centrohermitian perturbations that preserve the 
graph structure. We emphasize the important role of 
centrosymmetry or centrohermiticity here: it guarantees 
that the diagonal matrix elements have a constant value, 
which results in PST among the edge states of the net¬ 
work. 
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V 

Figure 3. Normalized distributions of the best efficiencies V of each realization of the the fermionic embedded Gaussian 

ensemble without centrosymmetry or centrohermiticity for l = 6: (a) fc = n = 1, (b) fc = l,n = 2, (c) k = 2, n = 2, 

(d) k = l,n = 3, (e) k = 2, nm 3, (f) k = n = 3, (g) k = l,n = 4, (h) k = 2,n = 4, (i) fc = 3,n = 4, (j) fc = n = 4 

(k) k = 1, n = 5, (1) k = 2, n = 5, (m) fc = 3, n = 5, (n) fc = 4, n = 5, (o) fc = n = 5. (Rows have the same particle number n, 

and columns the same fc value). Note that the two last rows correspond to filling factors above half-filling. Empty histograms 
correspond to /? = 1 and shaded (violet) to j3 = 2. 


IV. STATISTICS OF THE TRANSPORT 
EFFICIENCY FOR FERMIONS 


We consider now the case of the fermionic embedded 
Gaussian ensemble. For concreteness we analyze the case 
l = 6, varying the number of fermions n from 1 to 5, and 


rank of interaction k also from 1 to 5. In general, the 
distribution of efficiencies of the ensemble is rather broad 
with negligible contribution to efficiencies above 95% and 
thus, we will focus on distribution of best efficiencies V 
for each member of the ensemble. 

We begin describing first the case where no centrosym- 
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Figure 4. Same as Fig. [3] where we impose centrosymmetry or centrohermiticity in the one-particle space. Its extension to 
the k- and n-particle spaces is such that full centrosymmetry/centrohermiticity is achieved for odd values of k or n. Note 
that combining full centrosymmetry/centrohermiticity and being above half-filling, for k ~ 1/2 yields the best efficiencies per 
realization. 


metry or centrohermiticity is considered. The distribu¬ 
tion of the best efficiencies V of each realization are illus¬ 
trated in Fig [3] In this figure, all frames in a row have the 
same number of fermions n , and all frames in a column 
the same value of k. As before, the empty histograms 
correspond to (3 = 1 and the filled ones to /? = 2. In 
general, the lack of centrosymmetry or centrohermitic¬ 
ity yields poor efficiencies, with a marginal improvement 


for /3 = 2 compared to the (3 = 1 case. Interestingly, 
below and at half-filling (first three rows) only the case 
k = n = 1 marginally reaches the benchmark value, while 
above half-filling (last two rows of the figure) the distri¬ 
butions of V yield better results. In particular, the case 
n = 5 = l — 1 is the only one where some realizations 
achieve the benchmark value. Yet, the probability of 
such events is very small. 

































































The distributions of the best efficiencies of each real¬ 
ization including centrosymmetry/centrohermiticity are 
displayed in Fig. [4] As in the bosonic case, [3 = 1 seems 
to yield better efficiencies than the (3 = 2 case. For k = 1, 
the distributions of V for all n surpass the benchmark 
value, and have a non-zero probability of displaying PST, 
specially for odd fermion number. Yet, the k = 1 case 
does not correspond to the overall optimal situation: the 
probability that V is above the benchmark value is clearly 
larger for n = 5 = l — 1 and 2 < k < 4. Indeed, while 
we find that nearly 30% of the best efficiencies are larger 
than the benchmark for k = n = 1 and (3 = 1 (similar 
values are obtained for k = 1 and k = 5 for n = 5), for 
k = 3 and n = 5 this percentage is close to 59%. For 
(3 = 2 similar results are found with smaller percentages. 

While centrosymmetry clearly enhances the efficien¬ 
cies, its role is more involved. As we mentioned above, 
centrosymmetry/centrohermiticity is introduced at the 
one-particle level and then extended onto the k- and n- 
particle spaces using Eq. [I] The corresponding Jfc and 
J n matrices are exchange matrices or partial exchange 
matrices depending on the parameters k and n, and l. 
For the present case, l = 6, it can be shown that for odd 
particle-number space (k or n ) the corresponding ma¬ 
trix (Jfc or J n ) is a bona fide exchange matrix, whereas 
for even-values there are more than one states that are 
mapped onto itself, i.e., the corresponding J*, or J n ma¬ 
trix is a partial exchange matrix. Therefore, independent 
of fc, including the case where Jk is a partial exchange 
matrix, the n = 3 and n = 5 many-body Hamiltonians 
are fully centrosymmetric/centrohermitian. As shown 
in Fig. [4l full centrosymmetry or centrohermiticity, in 
comparison with the partial cases, yield the best effi¬ 
ciency distributions. When the many-body Hamiltonian 
is partial centrosymmetric/centrohermitian, the best dis¬ 
tributions correspond to the k = 1 case. In addition to 
centrosymmetry/centrohermiticity, the filling factor ap¬ 
pears to play an important role: above half-filling the 
distributions of V display the largest probability to yield 
efficiencies above the benchmark value; in terms of k , the 
optimal case corresponds to k ~ 1 /2. 

Regarding the pair of states that display the best ef¬ 
ficiencies when centrosymmetry/centrohermiticity is im¬ 
posed, as in the bosonic case, these occur among cen¬ 
trosymmetry related pairs, excluding 

those which are mapped onto themselves by J n when¬ 
ever J n is partial centrosymmetric. These pairs of states 
appear uniformly distributed, that is, there is no special 
pair of states that display better transport properties. 
This may be related to the fact that the fermionic graph 
is regular fli| . 

We have confirmed these results for l = 7, n = 6, k = 1, 
..., 6 and l = 8, n = 7, k = 1, ..., 7. 


V. CONCLUSIONS 

The primary aim of the present paper is to study 
transport of excitations in disordered networks with ran¬ 
dom fc-body interactions. This is important and cer¬ 
tainly of interest because of possible applications in quan¬ 
tum communication protocols [23| and artificial solar 
cells [24j]. Towards this end, we have studied the dis¬ 
tribution of quantum efficiencies in disordered networks 
with many-body interactions, whose structure is mod¬ 
eled by the embedded Gaussian ensemble, considering 
bosons and fermions, with and without time-reversal 
symmetry. In particular, we studied the role played by 
centrosymmetry/centrohermiticity, which is defined at 
the one-particle space, and then extended to the k- and 
n-particle spaces. We have shown that (full) centrosym¬ 
metry enhances the efficiencies dramatically, being a re¬ 
quirement to have non-zero probability for PST; the lack 
of centrosymmetry/centrohermiticity yields rather poor 
efficiencies. 

For bosons distributed in two single-particle levels, 
centrosymmetry is inherited in the k- and n-body spaces. 
In this case, PST is obtained for k = 1 for almost all real¬ 
izations; the fact that our computation of the efficiencies 
involves an upper bound for the time, T, constrains the 
relevant time scale for the achievement of the PST. In 
terms of fc, the probability of having PST decays with 
increasing k. However, we stress that for k > 1, the best 
efficient scenario is when k ~ nj 2, n is the total num¬ 
ber of bosons. With respect to the value of /3, the re¬ 
sults are marginally better when time-reversal symmetry 
is preserved. The pairs of states showing the best efficien¬ 
cies are those at the edges of the network, i.e., where all 
bosons are in one of the two single-particle levels. Then, 
in this case, state transfer corresponds to the physical 
transport of all bosons to the other single-particle state. 

For fermions, we found that full centrosymmetry/cen- 
trohermiticity of the n-particle Hamiltonian enhances 
considerably the best efficiencies, especially when the 
filling-factor is larger than 1/2. We note that cen- 
trosymmetry/centrohermiticity is indroduced at the one- 
particle level, and then extended to the n-body space. 
For fermions, the rank of the interaction which displays 
the highest probability that the best efficiency V is above 
the benchmark value corresponds to k ~ 1/2. A clear ex¬ 
planation of this is left open. The pairs of states that 
yield the best efficiencies appear uniformly among those 
states linked by centrosymmetry. 

Previous results have shown that random perturba¬ 
tions on networks with PST destroy or affect significantly 
this property; for details see (H, 26] and also [23] and 
references therein. Our results show that, despite of the 
random character of the fc-body interactions that we have 
considered, certain n-body networks display good effi¬ 
ciencies and may attain near perfect state transfer with 
non-zero probability. 

Our results could be exploited as new design princi¬ 
ples of networks with good efficiency, which is preserved 
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under certain many-body random perturbations. For in¬ 
stance, considering the implementation of efficient quan¬ 
tum wires, it may be interesting to consider the case of 
filling factors that are smaller but close to one, where 
many-body interactions yield robustly very good efficien¬ 
cies. Finally, the results in our paper open the possibility 
to understand the good efficiency properties experimen¬ 
tally observed in exciton transport in biological systems, 
such as the Fenna-Matthews-Olson complex [7|. |24| . 
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